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Abstract 

We study the macroscopic quantum tunneling of magnetization of the F = 1 spinor condensate 
interacting through dipole-dipole interaction with an external magnetic field applied along the lon- 
gitudinal or transverse direction. We show that the ground state energy and the effective magnetic 
moment of the system exhibit an interesting macroscopic quantum oscillation phenomenon origi- 
nating from the oscillating dependence of thermodynamic properties of the system on the vacuum 
angle. Tunneling between two degenerate minima are analyzed by means of an effective potential 
method and the periodic instanton method. 
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I. INTRODUCTION 

Rapid experimental progresses in realization of spinor condensates have generated 
fascinating opportunities to study the spin dynamics and magnetic properties of condensate 
atoms. The properties of spinor condensates under external magnetic field are investigated 
both experimentally and theoretically 4]. The spin-exchange interaction plays an im- 
portant role in these works, which is reminiscent of the exchange interaction responsible for 
interesting magnetic properties in solid. Other interactions, due to much weaker than the 
exchange interaction, are in most cases ignored. 

Since the spin degree of freedom becomes accessible in an optical trap, the magnetic 
dipole-dipole interaction which arises from the intrinsic or field-included magnetic dipole 
moment 0, |^ should be taken into account. The dipolar coupling was first considered 
between atoms at different sites, based on the assumption that the spin-exchange interaction 
dominates the on-site interaction [3]. Due to its long range and vectorial characters, more 
and more attentions are paid to this spinor dipolar condensate, for example, the ground 
state structure and spin dynamics were examined for this novel quantum system in a single 
trap Pi and in deep optical lattices ^j. 

Recently a major experimental breakthrough has been achieved in the condensation 
of chromium atoms ^^Cr which possesses a large magnetic dipole moment of QfiB {fJ'B is the 
Bohr magneton) in its ground state. The dipolar interaction in this condensate is a factor 
of 36 higher than that for alkali atoms, which makes possible the study of many dipolar 
phenomena and new kinds of quantum phase transitions predicted by theory. Indeed, the 
long-range and anisotropic magnetic dipole-dipole interaction in degenerate quantum gases 
has been shown to lead to an anisotropic deformation of the expanding Cr-BEC which 
depends on the orientation of the atomic dipole moments 

In the present paper we mainly consider a rich set of macroscopic quantum phenomena 
occurred in different quantum phases of spinor dipolar condensate. These phases may be 
tuned via modifying the trapping geometry so that various effective strengths of the dipolar 
interaction can be achieved. Based on the ground state structure of a spin-1 condensate 
with dipole-dipole interaction at zero temperature, macroscopic quantum tunneling and 
oscillations occur in different phases, which are quite similar to what happens in molecular 
magnets 
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The paper is organized as follows. In section 2, we introduce the dipolar spinor BEC 
model for T = under the single mode approximation, taking as a starting point the spin 
system analogous to the magnetization tunneling in magnetic particles. The ground state 
structure of the system is summarized in Section 3 for zero, longitudinal and transverse fields, 
respectively. Sections 4 and 5 are devoted to analyze the macroscopic quantum phenomena 
in the ground state for different phases under longitudinal or transversal external fields, in 
which cases tunneling between magnetizations arise naturally. Finally a brief summary is 
given in Section 6. 

II. SPINOR CONDENSATE WITH DIPOLAR INTERACTION 

Our starting point is the many-body Hamiltonian H proposed in Ref. f^, which describes 
a. F = 1 spinor condensate at zero temperature trapped in an axially symmetric potential 
Vext- Without loss of generality, the symmetry axis is conveniently chosen to be the quan- 
tization axis z. We consider here two atomic interaction terms, the short-range coUisional 
interaction and the long-range magnetic dipolar interaction, the competition of which gives 
rise to different quantum phases. Under an external magnetic field B, the second quantized 
Hamiltonian of the system reads 



where ■j/'a (r) (a = 0, ±1) are the field annihilation operators for an atom in the hyper- 
fine state \F = l,mp = a). The two coefficients Cq = 4-7?^^ (oq + 202) /3M and C2 = 
Anh"^ (a2 — clq) /3M characterize the density-density and spin-spin collisional interactions, 
respectively. Here a/ (/ = or 2) being the s-wave scattering length for spin-1 atoms 
in the combined symmetric channel of total spin /. The dipolar interaction parameter is 
Cd = yUo5'^/^B/4vr with g being Lande gf-factor. The Vext represent the external trapping po- 
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tential which is spin independent for a far off- resonant optical trap. And e = (r — r') /|r — r'| 
is a unit vector. In this study, the external field B is assumed to be spatially uniform. 

In order to simplify the Hamiltonian (Q), we usually adopt the single mode approximation 
(SMA): V'a(r) — 0(r)aQ, 14 1, where (r) is the spin independent spatial wave function 
of the condensate, is the annihilation operator for mp = a component. It is always 
safe to use this approximation for ferromagnetic interactions, however, it may break down 
for antiferromagnetic interactions if both the atomic number N and the magnetization A4 
are large. The interaction parameters in our case of study are such that |c2| <^ Cq and 
Cd ^ 0.1 |c2|. Under these conditions, the single mode approximation is expected to be 
valid if the trapping potential is axially symmetric. All integral terms in the long range 
interactions involving e^*'^'= vanish after integrating over the azimuthal angles of r — r'. We 
then obtain a much simpler Hamiltonian by assuming that the mode function (r) hence 
possesses the axial symmetry and dropping the spin- independent constant terms (tI: 

H = (4 - c',) V + 3c', [LI + no) - giiB^ ■ L, (2) 

where L = aJ^F^^/ja^ characterizes the total many-body angular momentum operator, 
^^0 = ^0^0 is the number operator for mp = ^ atoms. The new parameters are = 
(C2/2) / rfr|0 (r) |4 and c', = (q/4) / drdr'\(j) (r) (r') P (1 - 3 cos^ 6^) /|r- rf, with 6^ be- 
ing the polar angle of (r — r'). Eq. may be put into the following dimensionless form by 
rescaling it in energy unit Icg]: 

H = {±1- c) V + 3c (lI + no) - B' ■ L, (3) 



where -|- and — correspond to C2 > and C2 < 0, respectively. The new parameter c = 
c^/|c2| thus measures the relative strength of dipolar interaction with respect to the spin- 
exchange interaction. The dimensionless magnetic field B' = ^fyU^B/ |c2| = -B'(sin6', 0, cos^) 
is assumed to lie in the x2;-plane with an angle 9 relative to the z-axis. 



III. GROUND STATE STRUCTURE UNDER EXTERNAL FIELD 

We summarize the ground state structure which has been described in Ref. |7|]. In the 
absence of an external field, the ground state of our system is divided into three distinct 
regions A, B, and C in the C2-c^ parameter plane (see Figure 1). We denote here the 



4 



simultaneous eigenstate of and as |Z, m) with eigenvalues l{l-\-l)fi? and mh, respectively. 
In Region A (c^ > and c'^ > 4) the ground state is given by G — \N,{)), a quantum 
superposition of a chain of Fock states \N/2 — k,2k,N/2 — k) in which the numbers of 
atoms in the spins 1 and —1 are equal. In Region B (c'^ < and c'^ < — C2/2), G = \N, ±N) 
is a Fock state with all the population in either mj? = 1 or — 1. In these two regions, the no 
term is at least a factor oi 1/N smaller than the rest and therefore can be neglected safely. 
In Region C, however, the ground state is a little more complicated because the Uq term is 
expected to be important. In general it is expressed as G = ^i9i 1^? 0), a superposition of 
different angular momentum states with (^lJ^ = 0. 

When we apply a longitudinal field to the system along the 2;-axis (^^ = 0), the condensate 
in different parameter regions behave quite differently. In region A, the system can be 
mapped onto an easy-plane anisotropic particle with the transverse xy-plane being the easy- 
plane. The ground state is G — \N,m) with m depends on the field strength hnearly, with 
steps. Region B corresponds to a uniaxial anisotropic magnetic spin model with the easy 
axis z. The presence of the external magnetic field simply removes the two-fold degeneracy 
and forces the atoms into the fully polarized state G = \N, N). The ground state in Region 
G is changed into G = Xl/^m„ 9i K)''^o) while tuq increases with the field strength. 

Now we check the situation when a transverse field is applied along the x-axis, i.e., 
9 — 7r/2. Due to the easy plane anisotropy in Region A, the ground state of the condensate 
is fully polarized, however, in this case along the x-axis, G — \N, 0)^^^^. The situation in 
Region B is more interesting because the ground state here G = '^^gm \N,m) is two-fold 
degenerate and it provides another example that can exhibit macroscopic quantum tunneling 
of magnetization. Stepwise magnetization curve will appear in Region G - each step means 
the breaking of one spin singlet pair. 

For clarity, we choose two of above models which admit extensive study of tunneling 
of magnetization, i.e.. Region A with a longitudinal field and Region B with a transverse 
field. The dipolar spinor condensate thus provides another platform for the investigation of 
macroscopic quantum phenomena. 
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IV. MACROSCOPIC QUANTUM OSCILLATION OF MAGNETIZATION 



We first consider the spinor dipolar condensate under longitudinal field along 2;-axis in 
Region A, in which case the transverse xy-plane corresponds to an easy-plane anisotropy. 
After dropping the unimportant no and constant terms, the effective Hamiltonian can be 
obtained 

Hla = 3cLl - B'L,, (4) 

where c > 0. The model is precisely the same as that for a ferromagnetic particle 12, Q| 
with easy-plane anisotropy and a magnetic field along hard axis. The Hamiltonian is exactly 
diagonal in terms of the eigenstate and with the eigenvalue = 3cm^ — B'm. For zero 
temperature the magnetization increases stepwise as a consequence of the fact that m take 
only integer values. This model, on the other hand, provides a perfect manifestation of the 
B vacuum effect originating from the oscillating dependence of thermodynamic properties 
of the system on the "vacuum angle" The concept of the vacuum was developed 

mainly for the models of modern quantum field theory [l6| , but also in condensed media the 
nonperturbative vacuum is not a mathematical abstraction. The vacuum angle is nothing 
but the factor with which the total time derivative enters the Lagrangian. For Aharonov- 
! 3ohm problem in conductors with charge density waves, is the normalized magnetic flux 
IT^ . In the Josephson junction of mesoscopic sizes, the vacuum angle depends on the voltage 
applied to the junction We see shortly the magnetic field enters the Lagrangian and 
plays the role of vacuum angle in our spinor dipolar condensate system. 

We express the partition function as a spin coherent state path integral for large number 
of atoms N ^ 1 

Z = Trexp{-f3HLA) = J D {fi (n)} exp (Se) , (5) 

where the measure of the integration is decomposed into D{fi{n)} = 
Yl^^i {N sin 6 kd(j)kdOk/ 271). The semiclassical approximation of the partition function 
turns out to be the transition amplitude between two spin configurations conneted by 
periodic orbits with fixed imaginary time period (3. Following the usual procedure we 
represent the the state vector of the system as coherent states and slice the integral into Af 
identical pieces of length e = P/M . Inserting complete sets of states gives 

Z = {ufI exp {-[3Hla) |n/) = J] D {i^k (n)} J] (n^+il (1 - eHLA) \nk) (6) 

k=l k=0 



After some algebraic evaluation and finally passing to the time continuum limit M — > oo 
we obtain the Euclidean action in imaginary time r = it as (dots now denote r-derivatives) 

Se= dT[iN ^{1 -cos 9) + 3cN (N + 1)cos'^ 9- B'N cos 9] (7) 
Jo 

with f3 = l/ksT and T the temperature. Integrating over cos^ we map the magnetic system 
onto a particle problem with Lagrangian 

/: = ^ + ^e0, (8) 



where the effective mass ^e// = l/6c and = (1 — B' /6cN). This 



agrangian is exactly 



the one for O vacuum of the non-Abelian gauge field in references |16l Il8l|. We noticed that 
the second term of C is the total imaginary time derivative and has no effect on the classical 
equation of motion, while it indeed alerts the canonical momentum into II,^ = Tneff4> + ^0- 
In order to minimize the Euclidean action Se, we try to find the classical configuration, that 
is, the periodic instanton solution under the boundary condition 0„ (r + /5) = 0„ + 2nn. 
The result is 0„ = 2TmT/ (3 and the Euclidean action for this solution is Se = Sort^ + i2mTQ, 
where n is the winding number characterizing homotopically nonequivalent classes and Sq = 
Tc'^kBT/Sc. The Euclidean functional integral of the partition function contains thus an 
additional summation over the homotopic number and detailed calculation leads to 

oo 

Z= J2 ^n = ^3[e,exp(-so)], (9) 

n=— oo 

where i)^ {v, q) is the Jacobi theta function oscillating with 0. By means of the well known 
asymptotics of the Jacobi theta function, the ground state energy can be shown clearly 
oscillating with 0, i.e. the external magnetic field B 

E, = -ksTlnZ = + -^{{0}}^ (10) 

12c ^rrieff 

where {{x}} is the difference between x and its nearest integer. According to this, the 
external magnetic field induces quantum oscillations in the dipolar spinor condensate. The 
period of oscillation is shown to be 6B = Qc'^/gjiB (see Figure 2). 

The topological term in the Lagrangian (jSJ leads to the oscillation behavior of our sys- 
tem. We emphasize here the ground state (or in the language of quantum field theory, the 
instanton vacuum) of the condensate acquires a vacuum angle owing to the presence of the 
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external magnetic field, which breaks the symmetry and changes the topology of the sys- 
tem. In a charge density wave ring-shaped conductor placed in an external vector potential 
field the magnetic susceptibility and the electrical conductivity oscillate as a function of 
the flux while the voltage applied on the Josephson junction induces the oscillation 
of the effective capacitance of the junction In our system of investigation the macro- 
scopic quantum effect manifests itself by the oscillation of the effective magnetic momentum 
M = —ksTdln Z/dB as a function of the applied magnetic field. The magnetization curve 
at zero temperature is stepwise due to the macroscopic quantum oscillation of the effec- 
tive magnetic moment and the period of oscillation depends on the strength of the dipolar 
interaction uniquely. We show this oscillation period in Figure 2, together with the mathe- 
matical function {{x}}. The functional integral approach presented here becomes necessary 
when environmental friction or dissipation is included, which amounts to the introduction 
into the Euclidean Lagrangian of an additional term describing the non-local interaction of 
the instantons [2^. This is, however, beyond the scope of this paper. It is interesting to 
estimate the modulation period in the typical laboratory experiment on spinor condensates. 
For a condensate with ~ 5 x 10^ sodium atoms ^^Na and density p ~ 10^'^cm~^ [23] we 
obtain the regular interval of magnetic field 6B ~ lAGauss. This oscillation period can be 
as small as 1.5 x 10~^Gauss for 2 x 10^ rubidium atoms ^''Rb with density 2.6 x lO^^cm"^ as 
in the first experimental achievement of BEC in JILA [24] and can be as large as 10.5Gauss 
for 5 X 10^ chromium atoms ^^Cr with density lO^'^cm"^ as in the latest dipolar condensate 
experiments in Stuttgart HQ. Moreover, it can be adjusted flexibly by changing the 
trapping potential geometry. 



V. QUANTUM TUNNELING OF MAGNETIZATION BETWEEN TWO LOCAL 
MINIMA 

Now we prepare the system properly in Region B and apply a transverse field along the 
X-axis to the condensate. Again we drop the constant terms and are left with the following 
effective Hamiltonian 

Htb = -SdLl - B'L, (11) 

where d = \c\. The model describes a quantum spin system with the easy-axis anisotropy 
while the external field is along x-axis. This model has been extensively studied in the 
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context of spin tunneling. Classically under the influence of a weak transverse fleld along 
the X direction, the two energy minima move away from the zero flled positions {+z or 
—z) and towards x-axis while remaining in the xz plane. For < < QdN ^ they are 
located at 7_ = arcsin (i?'/6(iA^) and 7+ = tt — respectively, with the angle 7 between 
L and 2;-axis. The degeneracy is removed when B' > B'^^^ = 6dN, where the system is 
completely polarized by the external fleld and the two minima merge along x-axis. Quan- 
tum mechanically, the degeneracy is lifted before the magnetic fleld reaches B'^^^ due to the 
magnetization tunneling. A well-known consequence of the tunneling between two degener- 
ate states is the lifting of their degeneracy: The two new eigenstates are a symmetric and 
an antisymmetric superposition of the original states characterized by an energy difference 
(or tunneling splitting) AEq inversely proportional to the tunneling rate. The quantity of 
interest to determine the occurrence of tunneling is therefore this energy difference between 
the two lowest eigenstates of the Hamiltonian. 

To calculate analytically this energy splitting, we use the effective potential method 
which maps the spin system onto a particle system and the result can be easily obtained 
with the periodic instanton method. The Schrodinger equation Htb = I^P) in the Lz 
representation takes the form 

2{E + Srfm^) C„ + B'^/{N - m) {N + m + l)C™+i 

+B',/{N + m){N-m + l)Cm-i = (12) 

where m = — A^, — -|- 1, ■ ■ ■ A^ and Cm = for |m| > A^. 
Let us introduce the generating function, 

N ^ 

$= > — ™ = exp imx) . (13) 

±^^^{N-m)\{N + m)\ 



Multiplying equation (fT^ by a factor exp {mx) / [N — m)\ {N + m)\ and summing all 
terms with — A^ < m < N, we can transform the equation p2|l into 

3d<i>" - B' sinh x^' + {E + B'N cosh x) $ = 0. (14) 

In order to remove the flrst derivative term, let us deflne a new function \1/ = 
$ exp( — coshx). When we replace $ with the new Schrodinger equation can be 
written after dividing by A^^ in the form 

N'^^" + xij(K-U) = 0, (15) 



where the corresponding parameter k describe the dimensionless energy k = E/3dN'^, and 
U = (a cosh a; — 1)^ — is the effective potential well with a = B'/6dN. 

The value plays the role of the Planck constant h, the potential takes the form of 
a double well for a < 1, i.e., for external magnetic field not exceeding the saturation value 
^'sat ~ QdN. The two local minima thus play the role of the degenerate classical states and 
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the energy splitting of the lower states takes the following form according to 

AEn = AEog"/n! (16) 
with the splitting for the ground state 



AEo = ^huj/TTCe^^i-S/h). (17) 

Here Se is the Euclidean action evaluated along the trajectory from the left minimum 
X- = — cosh^^ i^/o) to the right x+ = cosh~^ (l/'^); ^ = vT^-o^ is the small oscillator 
frequency near the bottom of potential well x±. The asymptotic form of the instanton 
trajectory determined the constant C and q = C^/2Jiuj. For the potential U we have 



1 + VT 



SE = \n{ "^"'^ ) - VT^. (18) 



a 



Including the prefactor we have finally 



A^n = -^N^'^— ^ 



2\2N 



exp 



{2SeN + 1)V1^ . (19) 



Experimentally this level splitting can be measured by means of the resonance mea- 
surement developed by Awschalom j2^. Tunneling between two degenerate orientations of 
magnetization leads to the splitting of the "non-tunneling" ground state energy level into two 
levels separated by AEo. Correspondingly a very weak ac field of frequency AEq/H, will in- 
duce transitions between the two levels, which should result in the resonant absorption of the 
energy of the field. The atoms in the condensate are utterly identical and we do not have the 



(fTUj) is expressed 
the dipole-dipole 



problem of distribution of particle sizes and shapes. The level splittin g ec . 
in units of It is easily shown that for the sodium condensate in j22| 
interaction, i.e., the anisotropic energy in our model, is estimated as = 1.69 x 10^ Hz or 
11. 5/ii('^a quantity much smaller than the anisotropy energy of molecular magnets MriuAc 
or Fcg j26|. The level splitting can be greatly enhanced by a smaller number of atoms in 
the condensate and also by a stronger dipolar-dipolar interaction. Taking, as an example, 
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sodium atoms {N — 38) under an external field oi a — 0.6, we have AEq — 1.14 x lO'^Hz. 
For the condensate of ^^Cr, the anisotropic energy is enhanced to = 2.44 x lO^Hz or 
166/^i^. For N — 39, we have AE'q = 3.15 x 10^ Hz. The level splitting remains in the same 
magnitude order because it depends very sensitively on the total number of atoms. These 
data are easily accessible in the present ultracold atom experiments. 

VI. SUMMARY 

Inspired by the macroscopic quantum tunneling in the magnetic system, we have in- 
vestigated the macroscopic quantum tunneling in the dipolar spinor condensates at zero 
temperature and obtained some interesting results by analyzing different phase areas with 
applied external fields. We found that the ground state energy and the effective magnetic 
moment oscillate with the external magnetic field in Region A under a longitudinal field and 
the oscillating period depends on the strength of the dipolar interaction as SB — Qc'Jg^B- 
This model provides a condensed media realization of the © vacuum in quantum field theory. 
The model in Region B with a transverse field provides an example where quantum tunnel- 
ing of magnetization occurs between two local minima. We estimated the level splitting to 
be at the reach of current ultracold atom experiments. 
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Figure Captions: 

Figure 1: (Color online) Magnetic phase diagram of dipolar spinor condensate 
parametrized in the c'a-c^ plane. Corresponding ground states are shown for zero exter- 
nal field. The two tunneling models studied in this paper are in phase A with a longitudinal 
field and phase B with a transverse field. 

Figure 2: The mathematical function {{x}} and the oscillation of the magnetization with 
period SB. 
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